--------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\yiranli\Dropbox\A&A\Stata14\01_analysis.log
  log type:  text
 opened on:  11 Jan 2023, 10:24:52

. use 01_main.dta

. 
. 
. **============Table 1====================================
. sum econ misconduct lifebehavior behavior

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
        econ |        207           1           0          1          1
  misconduct |        112           1           0          1          1
lifebehavior |         43           1           0          1          1
    behavior |        118           1           0          1          1

. 
. **============Figure 1====================================
. tab chuliwithin if type==1

whether the |
        gov |
   responds |
 within one |
      month |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        134       65.37       65.37
          1 |         71       34.63      100.00
------------+-----------------------------------
      Total |        205      100.00

. tab chuliwithin if type==2

whether the |
        gov |
   responds |
 within one |
      month |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |         65       65.00       65.00
          1 |         35       35.00      100.00
------------+-----------------------------------
      Total |        100      100.00

. tab chuliwithin if type==3

whether the |
        gov |
   responds |
 within one |
      month |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |          6       21.43       21.43
          1 |         22       78.57      100.00
------------+-----------------------------------
      Total |         28      100.00

. tab chuliwithin if type==4

whether the |
        gov |
   responds |
 within one |
      month |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |         77       71.30       71.30
          1 |         31       28.70      100.00
------------+-----------------------------------
      Total |        108      100.00

. 
. 
. tab discipline if type==1

 punishment |      Freq.     Percent        Cum.
------------+-----------------------------------
         -1 |          1        0.50        0.50
          0 |         10        4.98        5.47
          2 |          1        0.50        5.97
          4 |          2        1.00        6.97
          5 |          4        1.99        8.96
          6 |          6        2.99       11.94
          7 |        112       55.72       67.66
          8 |         34       16.92       84.58
          9 |         25       12.44       97.01
         10 |          6        2.99      100.00
------------+-----------------------------------
      Total |        201      100.00

. tab discipline if type==2

 punishment |      Freq.     Percent        Cum.
------------+-----------------------------------
         -1 |          3        3.00        3.00
          0 |         21       21.00       24.00
          1 |         19       19.00       43.00
          2 |          1        1.00       44.00
          3 |          5        5.00       49.00
          4 |          4        4.00       53.00
          5 |         34       34.00       87.00
          6 |          8        8.00       95.00
          7 |          5        5.00      100.00
------------+-----------------------------------
      Total |        100      100.00

. tab discipline if type==3

 punishment |      Freq.     Percent        Cum.
------------+-----------------------------------
         -1 |          1        3.57        3.57
          0 |          2        7.14       10.71
          1 |          1        3.57       14.29
          2 |          1        3.57       17.86
          3 |          2        7.14       25.00
          4 |          2        7.14       32.14
          5 |          7       25.00       57.14
          6 |          2        7.14       64.29
          7 |          8       28.57       92.86
         10 |          2        7.14      100.00
------------+-----------------------------------
      Total |         28      100.00

. tab discipline if type==4

 punishment |      Freq.     Percent        Cum.
------------+-----------------------------------
         -1 |          1        0.93        0.93
          0 |         48       44.86       45.79
          1 |          5        4.67       50.47
          2 |          8        7.48       57.94
          3 |          2        1.87       59.81
          4 |          1        0.93       60.75
          5 |         21       19.63       80.37
          6 |          3        2.80       83.18
          7 |         18       16.82      100.00
------------+-----------------------------------
      Total |        107      100.00

. 
. 
. 
. **============Table 2====================================
. sum tt_repost_ln chuliwithin discipline type rank_code position_code

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
tt_repost_ln |        441    8.044858    3.360535          0   16.25869
 chuliwithin |        441    .3605442    .4807037          0          1
  discipline |        436    4.876147    3.135012         -1         10
        type |        441    2.088435    1.226183          1          4
   rank_code |        439    3.398633    2.006405          0          8
-------------+---------------------------------------------------------
position_c~e |        436     4.43578    1.668644          1          6

. sum econ misconduct lifebehavior behavior

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
        econ |        207           1           0          1          1
  misconduct |        112           1           0          1          1
lifebehavior |         43           1           0          1          1
    behavior |        118           1           0          1          1

. 
. tab rank_code,mi

  rank code |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |         40        9.07        9.07
          1 |         40        9.07       18.14
          2 |         82       18.59       36.73
          3 |         57       12.93       49.66
          4 |         88       19.95       69.61
          5 |         64       14.51       84.13
          6 |         39        8.84       92.97
          7 |         21        4.76       97.73
          8 |          8        1.81       99.55
          . |          2        0.45      100.00
------------+-----------------------------------
      Total |        441      100.00

. tab position_code,mi

   position |
       code |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |         17        3.85        3.85
          2 |         78       17.69       21.54
          3 |         43        9.75       31.29
          4 |         31        7.03       38.32
          5 |         94       21.32       59.64
          6 |        173       39.23       98.87
          . |          5        1.13      100.00
------------+-----------------------------------
      Total |        441      100.00

. tab chuliwithin

whether the |
        gov |
   responds |
 within one |
      month |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        282       63.95       63.95
          1 |        159       36.05      100.00
------------+-----------------------------------
      Total |        441      100.00

. tab discipline,mi

 punishment |      Freq.     Percent        Cum.
------------+-----------------------------------
         -1 |          6        1.36        1.36
          0 |         81       18.37       19.73
          1 |         25        5.67       25.40
          2 |         11        2.49       27.89
          3 |          9        2.04       29.93
          4 |          9        2.04       31.97
          5 |         66       14.97       46.94
          6 |         19        4.31       51.25
          7 |        143       32.43       83.67
          8 |         34        7.71       91.38
          9 |         25        5.67       97.05
         10 |          8        1.81       98.87
          . |          5        1.13      100.00
------------+-----------------------------------
      Total |        441      100.00

. 
. //Note: the detailed information on rainfall variable will be shown in 02_anal
> ysis.do
. 
. **============Table 3====================================
. probit chuliwithin tt_repost_ln i.type,vce(robust)

Iteration 0:   log pseudolikelihood = -288.29523  
Iteration 1:   log pseudolikelihood = -271.55923  
Iteration 2:   log pseudolikelihood = -271.53563  
Iteration 3:   log pseudolikelihood = -271.53563  

Probit regression                               Number of obs     =        441
                                                Wald chi2(4)      =      29.09
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -271.53563               Pseudo R2         =     0.0581

------------------------------------------------------------------------------
             |               Robust
 chuliwithin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
tt_repost_ln |   .0586075   .0194314     3.02   0.003     .0205226    .0966923
             |
        type |
          2  |   .0561453   .1595368     0.35   0.725    -.2565411    .3688317
          3  |   1.124811   .2847296     3.95   0.000      .566751    1.682871
          4  |  -.2041323   .1572051    -1.30   0.194    -.5122486     .103984
             |
       _cons |  -.8718917   .1829761    -4.77   0.000    -1.230518   -.5132651
------------------------------------------------------------------------------

. 
. probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code,vce(robust)

Iteration 0:   log pseudolikelihood =  -283.8899  
Iteration 1:   log pseudolikelihood = -262.35771  
Iteration 2:   log pseudolikelihood = -262.33873  
Iteration 3:   log pseudolikelihood = -262.33873  

Probit regression                               Number of obs     =        435
                                                Wald chi2(17)     =      39.83
                                                Prob > chi2       =     0.0014
Log pseudolikelihood = -262.33873               Pseudo R2         =     0.0759

-------------------------------------------------------------------------------
              |               Robust
  chuliwithin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0657331   .0209327     3.14   0.002     .0247058    .1067604
              |
         type |
           2  |   .0798487   .1691488     0.47   0.637     -.251677    .4113743
           3  |   1.103291    .295122     3.74   0.000     .5248626     1.68172
           4  |  -.2290298   .1724753    -1.33   0.184    -.5670751    .1090155
              |
position_code |
           2  |   .2960357   .4857153     0.61   0.542    -.6559488     1.24802
           3  |   .7271249   .4745105     1.53   0.125    -.2028986    1.657148
           4  |   .0182451   .4868807     0.04   0.970    -.9360237    .9725138
           5  |   .2619258   .4832585     0.54   0.588    -.6852434    1.209095
           6  |   .3698887   .4627121     0.80   0.424    -.5370102    1.276788
              |
    rank_code |
           1  |  -.2537306   .3648702    -0.70   0.487    -.9688631    .4614019
           2  |  -.0580815   .3385436    -0.17   0.864    -.7216148    .6054518
           3  |  -.0804885   .3508424    -0.23   0.819    -.7681269      .60715
           4  |  -.2958604   .3375749    -0.88   0.381    -.9574951    .3657743
           5  |  -.0577795   .3497844    -0.17   0.869    -.7433443    .6277852
           6  |  -.2906249    .377412    -0.77   0.441    -1.030339     .449089
           7  |  -.2581018   .4183191    -0.62   0.537    -1.077992    .5617885
           8  |   .1596154   .5983744     0.27   0.790    -1.013177    1.332408
              |
        _cons |  -1.120035   .3825346    -2.93   0.003    -1.869789   -.3702807
-------------------------------------------------------------------------------

. 
. oprobit discipline tt_repost_ln i.type,vce(robust)

Iteration 0:   log pseudolikelihood = -877.60027  
Iteration 1:   log pseudolikelihood = -749.66123  
Iteration 2:   log pseudolikelihood = -747.09507  
Iteration 3:   log pseudolikelihood = -747.08329  
Iteration 4:   log pseudolikelihood = -747.08329  

Ordered probit regression                       Number of obs     =        436
                                                Wald chi2(4)      =     165.39
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -747.08329               Pseudo R2         =     0.1487

------------------------------------------------------------------------------
             |               Robust
  discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
tt_repost_ln |   .0291965   .0160536     1.82   0.069    -.0022679    .0606609
             |
        type |
          2  |  -1.942705   .1659221   -11.71   0.000    -2.267906   -1.617503
          3  |  -1.218037    .291088    -4.18   0.000    -1.788559   -.6475145
          4  |  -2.076961   .1829882   -11.35   0.000    -2.435611   -1.718311
-------------+----------------------------------------------------------------
       /cut1 |  -3.590862   .3085869                     -4.195681   -2.986043
       /cut2 |  -2.062762   .2226773                     -2.499202   -1.626323
       /cut3 |  -1.828465   .2142148                     -2.248319   -1.408612
       /cut4 |  -1.727938   .2121938                      -2.14383   -1.312046
       /cut5 |  -1.646878   .2099885                     -2.058448   -1.235308
       /cut6 |  -1.567963   .2098027                     -1.979169   -1.156757
       /cut7 |  -.9665337   .1877142                     -1.334447   -.5986208
       /cut8 |  -.7712375   .1794406                     -1.122935   -.4195404
       /cut9 |   .7439869   .1577401                      .4348219    1.053152
      /cut10 |   1.267347   .1649449                      .9440605    1.590633
      /cut11 |   1.996692   .2172389                      1.570912    2.422473
------------------------------------------------------------------------------

. 
. oprobit discipline tt_repost_ln i.type i.position_code i.rank_code,vce(robust)

Iteration 0:   log pseudolikelihood = -862.14766  
Iteration 1:   log pseudolikelihood = -723.96975  
Iteration 2:   log pseudolikelihood = -721.07168  
Iteration 3:   log pseudolikelihood = -721.05672  
Iteration 4:   log pseudolikelihood = -721.05672  

Ordered probit regression                       Number of obs     =        430
                                                Wald chi2(17)     =     208.13
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -721.05672               Pseudo R2         =     0.1637

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0365783   .0171558     2.13   0.033     .0029535    .0702031
              |
         type |
           2  |  -1.945207   .1701366   -11.43   0.000    -2.278668   -1.611745
           3  |  -1.313129    .292167    -4.49   0.000    -1.885766   -.7404922
           4  |  -2.164254   .1946915   -11.12   0.000    -2.545842   -1.782666
              |
position_code |
           2  |   .7405002   .5201029     1.42   0.155    -.2788826    1.759883
           3  |    1.11508   .4976754     2.24   0.025      .139654    2.090506
           4  |   .7788076   .4964432     1.57   0.117    -.1942033    1.751818
           5  |   .6549976    .511233     1.28   0.200    -.3470006    1.656996
           6  |   1.020992   .5034802     2.03   0.043     .0341887    2.007795
              |
    rank_code |
           1  |   .1392049    .346756     0.40   0.688    -.5404244    .8188342
           2  |  -.1768154   .3301304    -0.54   0.592    -.8238591    .4702283
           3  |   -.128926   .3355767    -0.38   0.701    -.7866442    .5287922
           4  |  -.1546041   .3334537    -0.46   0.643    -.8081613    .4989532
           5  |  -.0858047   .3454609    -0.25   0.804    -.7628956    .5912862
           6  |  -.1095871    .362528    -0.30   0.762    -.8201288    .6009547
           7  |  -.1388389   .4192434    -0.33   0.741    -.9605409    .6828631
           8  |   .8545065   .4857698     1.76   0.079     -.097585    1.806598
--------------+----------------------------------------------------------------
        /cut1 |  -2.899588   .4579112                     -3.797077   -2.002098
        /cut2 |  -1.298896   .4332492                     -2.148049   -.4497436
        /cut3 |  -1.065305   .4314318                     -1.910895   -.2197139
        /cut4 |  -.9668406   .4314006                      -1.81237    -.121311
        /cut5 |  -.8806044   .4310521                     -1.725451   -.0357579
        /cut6 |  -.7977155   .4307142                       -1.6419    .0464687
        /cut7 |   -.172897   .4253664                       -1.0066    .6608059
        /cut8 |   .0179466   .4235682                     -.8122317     .848125
        /cut9 |   1.566588   .4167977                        .74968    2.383497
       /cut10 |   2.106161   .4218252                      1.279399    2.932923
       /cut11 |   2.870335    .449283                      1.989756    3.750913
-------------------------------------------------------------------------------

. 
. 
. 
. 
. **============Table 4====================================
. 
. oprobit discipline tt_repost_ln jine_ln if type==1,vce(robust)    

Iteration 0:   log pseudolikelihood =   -229.663  
Iteration 1:   log pseudolikelihood = -215.36806  
Iteration 2:   log pseudolikelihood = -215.29859  
Iteration 3:   log pseudolikelihood = -215.29855  
Iteration 4:   log pseudolikelihood = -215.29855  

Ordered probit regression                       Number of obs     =        164
                                                Wald chi2(2)      =      26.26
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -215.29855               Pseudo R2         =     0.0625

------------------------------------------------------------------------------
             |               Robust
  discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
tt_repost_ln |   .0780085   .0240039     3.25   0.001     .0309616    .1250553
     jine_ln |   .1522675   .0421449     3.61   0.000      .069665    .2348701
-------------+----------------------------------------------------------------
       /cut1 |  -1.102294   .4681186                      -2.01979   -.1847989
       /cut2 |  -.2309951    .342361                     -.9020104    .4400203
       /cut3 |  -.0932535   .3264491                     -.7330821     .546575
       /cut4 |   .0635624   .2996356                     -.5237127    .6508374
       /cut5 |   .2226296   .2938236                     -.3532541    .7985133
       /cut6 |    1.98964   .3162058                      1.369888    2.609392
       /cut7 |   2.700438   .3344664                      2.044896     3.35598
       /cut8 |   3.845676   .4153952                      3.031516    4.659836
------------------------------------------------------------------------------

. oprobit discipline tt_repost_ln jine_ln i.position_code i.rank_code if type==1
> ,vce(robust)

Iteration 0:   log pseudolikelihood = -229.02633  
Iteration 1:   log pseudolikelihood = -200.02839  
Iteration 2:   log pseudolikelihood = -199.71324  
Iteration 3:   log pseudolikelihood = -199.71294  
Iteration 4:   log pseudolikelihood = -199.71294  

Ordered probit regression                       Number of obs     =        163
                                                Wald chi2(15)     =      89.75
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -199.71294               Pseudo R2         =     0.1280

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |    .076174   .0266911     2.85   0.004     .0238604    .1284875
      jine_ln |   .1383697   .0437634     3.16   0.002     .0525951    .2241443
              |
position_code |
           2  |    .359424   .9254237     0.39   0.698    -1.454373    2.173221
           3  |   .3160951   .9202725     0.34   0.731    -1.487606    2.119796
           4  |   .4188177   .9295715     0.45   0.652    -1.403109    2.240744
           5  |    .156726   .9032423     0.17   0.862    -1.613596    1.927048
           6  |   .1490803    .867332     0.17   0.864    -1.550859     1.84902
              |
    rank_code |
           1  |   1.717037   .9326365     1.84   0.066    -.1108966    3.544971
           2  |   1.243215   .8668908     1.43   0.152    -.4558597     2.94229
           3  |    1.40322   .8797621     1.59   0.111    -.3210817    3.127522
           4  |   1.444632   .8807398     1.64   0.101    -.2815865     3.17085
           5  |    1.37447   .8609306     1.60   0.110    -.3129235    3.061863
           6  |   1.542264   .8796979     1.75   0.080    -.1819124     3.26644
           7  |   1.639735   1.013364     1.62   0.106    -.3464215    3.625891
           8  |   2.917711   1.010081     2.89   0.004     .9379883    4.897433
--------------+----------------------------------------------------------------
        /cut1 |  -.0117417   .5033519                     -.9982932    .9748098
        /cut2 |   1.019732   .4235307                      .1896273    1.849837
        /cut3 |   1.214492   .4351652                      .3615843      2.0674
        /cut4 |   1.420201   .4238848                      .5894021       2.251
        /cut5 |   1.606422   .4342962                      .7552171    2.457627
        /cut6 |   3.525656   .4789266                      2.586977    4.464335
        /cut7 |   4.313286    .493434                      3.346173    5.280399
        /cut8 |   5.584229   .5406589                      4.524557    6.643901
-------------------------------------------------------------------------------

.           
. 
. char   xingchouwen[omit] 2 

. xi:oprobit discipline tt_repost_ln i.xingchouwen  if lifebehavior==1,vce(robus
> t) 
i.xingchouwen     _Ixingchouw_1-3     (naturally coded; _Ixingchouw_2 omitted)

Iteration 0:   log pseudolikelihood = -73.589349  
Iteration 1:   log pseudolikelihood = -68.772069  
Iteration 2:   log pseudolikelihood = -68.759809  
Iteration 3:   log pseudolikelihood = -68.759806  

Ordered probit regression                       Number of obs     =         40
                                                Wald chi2(3)      =       6.61
                                                Prob > chi2       =     0.0855
Log pseudolikelihood = -68.759806               Pseudo R2         =     0.0656

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .1567036   .1037493     1.51   0.131    -.0466413    .3600485
_Ixingchouw_1 |   .3022781   .5834656     0.52   0.604    -.8412935     1.44585
_Ixingchouw_3 |   .8454181   .5917361     1.43   0.153    -.3143635      2.0052
--------------+----------------------------------------------------------------
        /cut1 |   .0551151   1.187861                     -2.273049    2.383279
        /cut2 |   .2980116   1.211542                     -2.076568    2.672591
        /cut3 |   .4988106   1.242587                     -1.936615    2.934236
        /cut4 |   .8040762   1.284825                     -1.714134    3.322287
        /cut5 |   1.034574   1.288191                     -1.490233    3.559382
        /cut6 |   1.655098    1.34107                     -.9733501    4.283546
        /cut7 |   1.808188   1.357053                     -.8515865    4.467963
        /cut8 |   3.313774    1.54843                      .2789076    6.348641
        /cut9 |   3.623626   1.606913                      .4741338    6.773118
       /cut10 |   3.824472   1.694616                      .5030853    7.145859
-------------------------------------------------------------------------------

. xi:oprobit discipline tt_repost_ln i.xingchouwen i.position_code i.rank_code i
> f lifebehavior==1,vce(robust)
i.xingchouwen     _Ixingchouw_1-3     (naturally coded; _Ixingchouw_2 omitted)
i.position_code   _Iposition__1-6     (naturally coded; _Iposition__1 omitted)
i.rank_code       _Irank_code_0-8     (naturally coded; _Irank_code_0 omitted)

note: _Irank_code_8 omitted because of collinearity
Iteration 0:   log pseudolikelihood = -73.589349  
Iteration 1:   log pseudolikelihood = -56.270851  
Iteration 2:   log pseudolikelihood = -55.603295  
Iteration 3:   log pseudolikelihood = -55.575499  
Iteration 4:   log pseudolikelihood = -55.572112  
Iteration 5:   log pseudolikelihood = -55.571584  
Iteration 6:   log pseudolikelihood = -55.571508  
Iteration 7:   log pseudolikelihood = -55.571493  
Iteration 8:   log pseudolikelihood =  -55.57149  

Ordered probit regression                       Number of obs     =         40
                                                Wald chi2(13)     =          .
                                                Prob > chi2       =          .
Log pseudolikelihood =  -55.57149               Pseudo R2         =     0.2448

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .3504767   .1077257     3.25   0.001     .1393383    .5616152
_Ixingchouw_1 |   1.536465   .6338845     2.42   0.015     .2940738    2.778855
_Ixingchouw_3 |   1.898096   .7065428     2.69   0.007     .5132978    3.282895
_Iposition__2 |   10.43118   1.032618    10.10   0.000     8.407287    12.45508
_Iposition__3 |   10.95139   1.446727     7.57   0.000     8.115861    13.78693
_Iposition__4 |   11.72586   1.594992     7.35   0.000     8.599731    14.85199
_Iposition__5 |   10.80343   1.332411     8.11   0.000     8.191948     13.4149
_Iposition__6 |   10.15986   .8529218    11.91   0.000     8.488162    11.83155
_Irank_code_1 |  -10.89866   1.121341    -9.72   0.000    -13.09645   -8.700875
_Irank_code_2 |  -10.46459   1.127291    -9.28   0.000    -12.67404   -8.255138
_Irank_code_3 |  -12.42572   1.532784    -8.11   0.000    -15.42992    -9.42152
_Irank_code_4 |  -10.60184    .896887   -11.82   0.000     -12.3597   -8.843972
_Irank_code_5 |  -11.17515   1.002981   -11.14   0.000    -13.14095   -9.209342
_Irank_code_6 |  -10.18226   1.239561    -8.21   0.000    -12.61175   -7.752761
_Irank_code_7 |  -12.01429   1.340916    -8.96   0.000    -14.64244   -9.386143
_Irank_code_8 |          0  (omitted)
--------------+----------------------------------------------------------------
        /cut1 |   1.866553    1.09963                     -.2886809    4.021788
        /cut2 |   2.194521   1.168202                     -.0951134    4.484156
        /cut3 |   2.467144   1.150748                      .2117199    4.722568
        /cut4 |   2.886696   1.125392                      .6809687    5.092424
        /cut5 |   3.179264   1.113577                      .9966939    5.361834
        /cut6 |   3.978724   1.102396                      1.818068     6.13938
        /cut7 |   4.176936   1.127272                      1.967524    6.386348
        /cut8 |   6.342579   1.322685                      3.750165    8.934993
        /cut9 |   6.903499   1.257276                      4.439284    9.367715
       /cut10 |   7.406819   1.331346                      4.797429    10.01621
-------------------------------------------------------------------------------
Note: 1 observation completely determined.  Standard errors questionable.

. 
. 
. **============Figure 2 [panel A]===========================
. 
. probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code if rank_cod
> e<=5

Iteration 0:   log likelihood = -240.75549  
Iteration 1:   log likelihood = -221.71174  
Iteration 2:   log likelihood = -221.69831  
Iteration 3:   log likelihood = -221.69831  

Probit regression                               Number of obs     =        368
                                                LR chi2(14)       =      38.11
                                                Prob > chi2       =     0.0005
Log likelihood = -221.69831                     Pseudo R2         =     0.0792

-------------------------------------------------------------------------------
  chuliwithin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0451899   .0220723     2.05   0.041     .0019291    .0884508
              |
         type |
           2  |   .1989246    .182702     1.09   0.276    -.1591647    .5570139
           3  |   1.226743   .2978273     4.12   0.000     .6430125    1.810474
           4  |  -.0336489   .1858252    -0.18   0.856    -.3978597    .3305618
              |
position_code |
           2  |   .2065563   .5060393     0.41   0.683    -.7852626    1.198375
           3  |   .6965678   .4917742     1.42   0.157    -.2672919    1.660428
           4  |   -.217561   .5217183    -0.42   0.677     -1.24011    .8049881
           5  |   .1008482   .5022517     0.20   0.841     -.883547    1.085243
           6  |   .2988928   .4759346     0.63   0.530    -.6339218    1.231707
              |
    rank_code |
           1  |  -.2347863   .3937094    -0.60   0.551    -1.006443    .5368699
           2  |     -.0307   .3606788    -0.09   0.932    -.7376175    .6762176
           3  |   .0138991   .3709024     0.04   0.970    -.7130562    .7408544
           4  |  -.2201501   .3698487    -0.60   0.552    -.9450404    .5047401
           5  |   .0915509   .3804879     0.24   0.810    -.6541917    .8372936
              |
        _cons |  -1.021028    .379994    -2.69   0.007    -1.765803   -.2762534
-------------------------------------------------------------------------------

. 
. oprobit discipline tt_repost_ln i.type i.position_code i.rank_code if rank_cod
> e<=5

Iteration 0:   log likelihood =  -703.8302  
Iteration 1:   log likelihood = -596.22717  
Iteration 2:   log likelihood = -594.06757  
Iteration 3:   log likelihood = -594.05777  
Iteration 4:   log likelihood = -594.05777  

Ordered probit regression                       Number of obs     =        363
                                                LR chi2(14)       =     219.54
                                                Prob > chi2       =     0.0000
Log likelihood = -594.05777                     Pseudo R2         =     0.1560

-------------------------------------------------------------------------------
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0357848   .0182884     1.96   0.050    -.0000598    .0716293
              |
         type |
           2  |  -1.880536   .1736968   -10.83   0.000    -2.220975   -1.540096
           3  |  -1.295963   .2338072    -5.54   0.000    -1.754217   -.8377096
           4  |  -2.040965   .1766033   -11.56   0.000    -2.387102   -1.694829
              |
position_code |
           2  |   .7570538     .40439     1.87   0.061     -.035536    1.549644
           3  |   1.127826   .3971074     2.84   0.005     .3495098    1.906142
           4  |   .4454608   .4127979     1.08   0.281    -.3636083     1.25453
           5  |   .6017769   .4025027     1.50   0.135     -.187114    1.390668
           6  |   1.083013   .3818709     2.84   0.005     .3345602    1.831467
              |
    rank_code |
           1  |   .0913295   .3099222     0.29   0.768    -.5161069    .6987659
           2  |  -.2229118   .2876698    -0.77   0.438    -.7867342    .3409106
           3  |  -.0734607    .298056    -0.25   0.805    -.6576397    .5107183
           4  |  -.1485129   .2907598    -0.51   0.610    -.7183916    .4213658
           5  |  -.0001928   .3023364    -0.00   0.999    -.5927613    .5923757
--------------+----------------------------------------------------------------
        /cut1 |  -2.887808   .3524318                     -3.578562   -2.197055
        /cut2 |  -1.295549   .3152116                     -1.913353   -.6777462
        /cut3 |  -1.039643   .3155599                     -1.658129   -.4211572
        /cut4 |  -.9284259   .3157436                     -1.547272   -.3095798
        /cut5 |   -.863562   .3158237                     -1.482565   -.2445589
        /cut6 |  -.8012795   .3158398                     -1.420314   -.1822449
        /cut7 |  -.1342425   .3127538                     -.7472287    .4787437
        /cut8 |   .0472219   .3108888                      -.562109    .6565527
        /cut9 |   1.753856   .3189042                      1.128815    2.378897
       /cut10 |   2.259932   .3295108                      1.614102    2.905761
       /cut11 |   2.823671   .3548653                      2.128148    3.519195
-------------------------------------------------------------------------------

. 
. 
. 
. **============Figure 2 [panel B]===========================
. 
. 
. probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code if  positio
> n_code==5 | position_code==6 | position_code==1 

Iteration 0:   log likelihood = -184.40508  
Iteration 1:   log likelihood = -167.93091  
Iteration 2:   log likelihood = -167.90841  
Iteration 3:   log likelihood = -167.90841  

Probit regression                               Number of obs     =        283
                                                LR chi2(14)       =      32.99
                                                Prob > chi2       =     0.0029
Log likelihood = -167.90841                     Pseudo R2         =     0.0895

-------------------------------------------------------------------------------
  chuliwithin |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0659311   .0257951     2.56   0.011     .0153735    .1164886
              |
         type |
           2  |   .0076727   .2174968     0.04   0.972    -.4186131    .4339586
           3  |   1.128288   .3578334     3.15   0.002     .4269477    1.829629
           4  |  -.2371198   .2172585    -1.09   0.275    -.6629387    .1886992
              |
position_code |
           5  |   .2490682   .6141309     0.41   0.685    -.9546063    1.452743
           6  |   .2917663   .5790821     0.50   0.614    -.8432138    1.426746
              |
    rank_code |
           1  |   .0235896   .5489905     0.04   0.966    -1.052412    1.099591
           2  |   .1679626   .4945899     0.34   0.734    -.8014157    1.137341
           3  |  -.1100196   .5354309    -0.21   0.837    -1.159445    .9394058
           4  |  -.5832184   .5242849    -1.11   0.266    -1.610798    .4443611
           5  |   .0137725   .5303168     0.03   0.979    -1.025629    1.053174
           6  |  -.2806677   .5598657    -0.50   0.616    -1.377984    .8166489
           7  |   .0472035   .6369261     0.07   0.941    -1.201149    1.295556
           8  |  -.2311194   .9002001    -0.26   0.797    -1.995479     1.53324
              |
        _cons |  -1.109566     .40131    -2.76   0.006    -1.896119   -.3230128
-------------------------------------------------------------------------------

. 
. oprobit discipline tt_repost_ln i.type i.position_code i.rank_code if  positio
> n_code==5 | position_code==6| position_code==1 

Iteration 0:   log likelihood = -554.54136  
Iteration 1:   log likelihood = -461.85974  
Iteration 2:   log likelihood = -459.85452  
Iteration 3:   log likelihood = -459.84561  
Iteration 4:   log likelihood = -459.84561  

Ordered probit regression                       Number of obs     =        282
                                                LR chi2(14)       =     189.39
                                                Prob > chi2       =     0.0000
Log likelihood = -459.84561                     Pseudo R2         =     0.1708

-------------------------------------------------------------------------------
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0511012   .0205283     2.49   0.013     .0108664    .0913361
              |
         type |
           2  |  -1.836713   .1983523    -9.26   0.000    -2.225476    -1.44795
           3  |  -1.226511   .2804499    -4.37   0.000    -1.776183   -.6768395
           4  |  -2.210695   .2043949   -10.82   0.000    -2.611301   -1.810088
              |
position_code |
           5  |   .8016493   .4855081     1.65   0.099    -.1499291    1.753228
           6  |   1.204129   .4594229     2.62   0.009     .3036766    2.104581
              |
    rank_code |
           1  |   .1180666   .4308826     0.27   0.784    -.7264478     .962581
           2  |  -.3461362   .3901328    -0.89   0.375    -1.110782    .4185099
           3  |  -.2157718   .4238398    -0.51   0.611    -1.046483    .6149389
           4  |  -.3825197   .4051744    -0.94   0.345    -1.176647    .4116076
           5  |  -.0778145   .4208083    -0.18   0.853    -.9025838    .7469547
           6  |  -.2263646   .4384052    -0.52   0.606    -1.085623    .6328937
           7  |  -.6693247   .5061241    -1.32   0.186     -1.66131    .3226603
           8  |    .644505   .7378558     0.87   0.382    -.8016659    2.090676
--------------+----------------------------------------------------------------
        /cut1 |  -2.791237   .3803706                     -3.536749   -2.045724
        /cut2 |  -1.262771   .3271188                     -1.903912   -.6216295
        /cut3 |  -.9570355   .3280608                     -1.600023   -.3140481
        /cut4 |  -.8339089   .3284155                     -1.477591   -.1902264
        /cut5 |  -.7726139   .3285131                     -1.416488   -.1287401
        /cut6 |  -.6846091   .3285216                       -1.3285   -.0407186
        /cut7 |   .0224757   .3243517                     -.6132421    .6581934
        /cut8 |    .185964   .3221365                     -.4454118    .8173399
        /cut9 |   1.896227   .3329844                       1.24359    2.548865
       /cut10 |     2.4089    .346332                      1.730102    3.087699
       /cut11 |    3.00278   .3768233                       2.26422     3.74134
-------------------------------------------------------------------------------

. 
. 
. 
. **============Table 5===========================
. 
. tobit discipline tt_repost_ln i.type i.position_code i.rank_code, ll(0) vce(ro
> bust) 

Tobit regression                                Number of obs     =        430
                                                F(  17,    413)   =      19.28
                                                Prob > F          =     0.0000
Log pseudolikelihood = -923.22635               Pseudo R2         =     0.1163

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .1151772   .0447928     2.57   0.010     .0271268    .2032275
              |
         type |
           2  |   -4.09122   .3517558   -11.63   0.000    -4.782675   -3.399765
           3  |  -2.524943   .5657866    -4.46   0.000    -3.637123   -1.412762
           4  |  -5.353156   .4469615   -11.98   0.000    -6.231759   -4.474552
              |
position_code |
           2  |   1.886958   1.518497     1.24   0.215    -1.097989    4.871904
           3  |   3.326023   1.486619     2.24   0.026      .403739    6.248306
           4  |   2.024828   1.551957     1.30   0.193    -1.025891    5.075547
           5  |   1.889885   1.496857     1.26   0.207    -1.052523    4.832294
           6  |   3.039978    1.45863     2.08   0.038     .1727131    5.907242
              |
    rank_code |
           1  |   .5022044   .8948203     0.56   0.575    -1.256766    2.261175
           2  |  -.5012383   .8274769    -0.61   0.545     -2.12783    1.125353
           3  |  -.3953237   .8961887    -0.44   0.659    -2.156984    1.366336
           4  |   -.534904    .854231    -0.63   0.532    -2.214087    1.144279
           5  |  -.2698759   .8766253    -0.31   0.758     -1.99308    1.453328
           6  |  -.8354083   .9107098    -0.92   0.360    -2.625613    .9547964
           7  |  -.5670609   1.024815    -0.55   0.580    -2.581566    1.447444
           8  |   .6567715   1.074228     0.61   0.541    -1.454865    2.768408
              |
        _cons |   3.988043   1.297481     3.07   0.002     1.437552    6.538534
--------------+----------------------------------------------------------------
       /sigma |   2.745189   .1333165                      2.483125    3.007252
-------------------------------------------------------------------------------
            86  left-censored observations at discipline <= 0
           344     uncensored observations
             0 right-censored observations

. 
. drop if discipline==0
(81 observations deleted)

. oprobit discipline tt_repost_ln i.type i.position_code i.rank_code,vce(robust)

Iteration 0:   log pseudolikelihood = -655.55875  
Iteration 1:   log pseudolikelihood = -535.11593  
Iteration 2:   log pseudolikelihood = -532.03744  
Iteration 3:   log pseudolikelihood = -532.02396  
Iteration 4:   log pseudolikelihood = -532.02396  

Ordered probit regression                       Number of obs     =        350
                                                Wald chi2(17)     =     201.93
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -532.02396               Pseudo R2         =     0.1884

-------------------------------------------------------------------------------
              |               Robust
   discipline |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
 tt_repost_ln |   .0431693   .0191827     2.25   0.024     .0055719    .0807668
              |
         type |
           2  |  -2.355766   .2057721   -11.45   0.000    -2.759072    -1.95246
           3  |  -1.567254   .3511123    -4.46   0.000    -2.255422   -.8790869
           4  |  -1.945249   .2181043    -8.92   0.000    -2.372725   -1.517772
              |
position_code |
           2  |   .0759679   .6060192     0.13   0.900    -1.111808    1.263744
           3  |   .1876158   .5662587     0.33   0.740     -.922231    1.297463
           4  |   .0633133   .6129547     0.10   0.918    -1.138056    1.264682
           5  |  -.1894382   .6050974    -0.31   0.754    -1.375407    .9965308
           6  |   .0586527   .5940816     0.10   0.921    -1.105726    1.223031
              |
    rank_code |
           1  |  -.1347412   .4927694    -0.27   0.785    -1.100551    .8310691
           2  |  -.3285878   .4753391    -0.69   0.489    -1.260235    .6030597
           3  |  -.1454598   .4824038    -0.30   0.763    -1.090954    .8000343
           4  |  -.2006175   .4799795    -0.42   0.676     -1.14136    .7401251
           5  |   -.209526   .4989695    -0.42   0.675    -1.187488    .7684362
           6  |  -.0675944    .515931    -0.13   0.896    -1.078801    .9436118
           7  |  -.1186035   .5581175    -0.21   0.832    -1.212494    .9752868
           8  |   .8761347   .6425401     1.36   0.173    -.3832207     2.13549
--------------+----------------------------------------------------------------
        /cut1 |  -3.760222   .4887566                     -4.718167   -2.802277
        /cut2 |  -2.868639   .4397997                     -3.730631   -2.006648
        /cut3 |  -2.638628   .4395389                     -3.500109   -1.777148
        /cut4 |  -2.468424   .4365897                     -3.324124   -1.612723
        /cut5 |  -2.322038   .4356537                     -3.175904   -1.468173
        /cut6 |  -1.391432   .4158205                     -2.206425   -.5764389
        /cut7 |  -1.132174   .4102426                     -1.936234   -.3281129
        /cut8 |   .6486478   .3947528                     -.1250534    1.422349
        /cut9 |   1.216437   .3995488                      .4333354    1.999538
       /cut10 |   1.985198   .4266125                      1.149053    2.821343
-------------------------------------------------------------------------------

. 
. 
. log close
      name:  <unnamed>
       log:  C:\Users\yiranli\Dropbox\A&A\Stata14\01_analysis.log
  log type:  text
 closed on:  11 Jan 2023, 10:24:55
--------------------------------------------------------------------------------
